Surface plasmon Fourier optics 
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Surface plasmons are usually described as surface waves with either a complex wavevector or a 
complex frequency. When discussing their merits in terms of field confinment or enhancement of the 
local density of states, controversies regularly arise as the results depend on the choice of a complex 
wavevector or a complex frequency. In particular, the shape of the dispersion curves depends on 
this choice. When discussing diffraction of surface plasmon a scalar approximation is often used. 
In this work, we derive two equivalent vectorial representations of a surface plasmon field using 
an expansion over surface waves with either a complex wavevector or a complex frequency. These 
representations can be used to account for propagation and diffraction of surface waves. They can 
also be used to discuss the issue of field confinment and local density of states as they have a 
non-ambiguous relation with the two dispersion relations. 
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I. INTRODUCTION 

Surface plasmons have been known since the pioneer- 
ing work of Ritchie in the 1950si. Considerable advances 
made in nanotechnology in recent years and the desire 
to control and manipulate light at nanoscale have 
renewed the interest in surface plasmons 2 -. Numerical 
simulations and experiments have demonstrated unique 
properties of different plasmonic nanostructure such 
as extraordinary transmission 3 -" 4 -, guiding 5 -* 6 -! 7 -^, fluores- 
cence enhancement ^ iV 1 ' 12 ' 13 , field enhancemen t 14 ' 15 ' 16 , 
focussing-i 7 -, superresolutio n 20 i 21 i 22 , omnidirectional 
absorptio n 23 ' 25 , coherent thermal emissio n 24 ! 25 ' 26 . 

In this paper, we shall focus on surface plasmons prop- 
agating along flat surfaces. Propagation of surface plas- 
mons on a flat surface perpendicular to the z axis is of- 
ten discussed using a mode E(z) exp[i(K x x + K y y — ut)] 
characterized by a frequency lu and a wave vector K = 
Kjjc + K y y parallel to the interface. However, the sur- 
face plasmon fields diffracted by edges, guided by ridges, 
focussed by lenses cannot be described by a simple mode. 
It is well-known that a finite size beam propagating in a 
vacuum has to be described in terms of a linear superpo- 
sition of plane waves. Differents ansatz, often neglecting 
polarization, have been used in the litterature to address 
this question 5 - ' 18 ' 27 . One of the goals of this paper is to 
derive a rigorous representation for the surface plasmon 
field. Such a superposition is the equivivalent of the an- 
gular plane wave spectrum for surface plasmons It can be 
used to develope a framework for surface plasmon fourier 
optics. 

In doing so, a difficulty arises. When losses are taken 
into account, a mode with real K and real uj is no longer 
a valid solution. Although we can still use a Fourier rep- 
resentation with real K and real ui, it is not convenient 
to deal with waves that are not a solution. Elementary 
solutions using either a complex K or a complex u can 
be found. However, we cannot assume that they form a 



basis. The first issue is thus to derive a general repre- 
sentation for the surface plasmon field as a superposition 
of modes. The second issue is related to the dispersion 
relation. A dispersion relation can be found when using 
either a complex K and a real u> or vice versa. These two 
choices leads to different shapes as seen in FigJTJ One 
dispersion relation has an asymptote for very large val- 
ues of K while the other has limited values of K and 
presents a backbending. 

This issue was first noted by Alexander et al^ 8 - and 
later discussed by Arakawa et ali^. Arakawa remarked 
that when plotting the position of the dips in a reflec- 
tivity experiment where the angle of incidence is varied 
at fixed frequency, one finds the dispersion relation with 
backbending. Instead, when plotting the points obtained 
from a spectrum at fixed angle, one finds the dispersion 
relation without backbending. This approach seems sim- 
ple and easily applicable. It is sufficient to explain the 
attenuated total reflection (ATR) experiments. Never- 
theless, it is not a general prescription that can be used 
to discuss all possible issues. Let us illustrate this point 
by addressing two questions regarding the most impor- 
tant properties of surface plasmons : confinment of the 
fields and large density of states. For a theoretical dis- 
cussion of these applications, different dispersion rela- 
tions lead to different predictions. Confinment of the 
field is the key property regarding applications such as 
optical lithography, enhanced non-linear effects or super- 
resolution issues. The dispersion relation with a back- 
bending predicts a cut-off spatial frequency and therefore 
a resolution limit whereas the dispersion relation with- 
out backbending does not predict any resolution limit. 
Enhancement of the local density of states (or Purcell 
effect) is fundamental for fluorescence enhancement and 
more light emission assisted by surface plasmons. The 
dispersion relation with a backbending again predicts a 
cut-off spatial frequency and therefore an upper limit to 
the LDOS. No limit is predicted by the other dispersion 
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FIG. 1: (color online) Dispersion of surface plasmons prop- 
agating along the metal/air interface, (a) Real to is chosen 
to obtain a complex K ap . (b) Real K is chosen to obtain a 
complex ui ap . The flat asymptote (dashed line) situated at 
u! p /\/2 represents the nonretarded surface plasmon solution. 
The slanting solid line represents the light cone inside which a 
wave is propagating (radiative) and outside which is evanes- 
cent (surface). 



relation. A general discussion on the applicability of the 
different dispersion relations is thus needed. 

In this paper we start by deriving two general represen- 
tations of the surface plasmon field in terms of linear su- 
perpositions of modes having the structure exp[i(K x x + 
K y y + 72: — cot)] with a well-defined polarization. Each 
representation is associated with either a complex fre- 
quency or a complex wavevector and therefore to a par- 
ticular dispersion relation. We then show that the most 
convenient choice depends on the physical problem to be 
discussed. We introduce a prescription that allows to 
choose complex or real frequency and the associated dis- 
persion relation. We then apply our analysis to discuss 
the resolution limit issue and the LDOS issue. The paper 



is organized as follows. For the sake of completeness, we 
briefly summarize the derivation of the dispersion rela- 
tion in the next section. The following section introduces 
the general representations of the surface plasmon field. 
We then discuss the physical content of these represen- 
tations. 



II. SURFACE PLASMON DISPERSION 
RELATION 

Let us consider a flat metal surface z = bounded by 
dielectric media with dielectric constant e\. For conve- 
nience, we describe the dielectric response of the metal 
to an electric field using the local Drude model 



e 2 (w) = 1 - 



lu{lo + iv e ) ' 



(1) 



where lo p is the bulk plasmon frequency and v e is a phe- 
nomenological bulk electron relaxation rate. We derive 
the dispersion relation of surface plasmons propagating 
along the metal/dielectric interface. 

We search a solution of Maxwell equations for an inter- 
face between two linear isotropic and local media charac- 
terized by dielectric constants e TO where m = 1, 2 denotes 
medium 1 [z < 0) or 2 (z > 0). A surface wave solution 
has a structure exp[z(K ■ r + j m \z\ — uit)] with 
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where 7 m is chosen so that 3m7 m > 0. Boundary condi- 
tions impose the continuity of the tangential components 
of the electric field and eE z . It follows that a p-polarized 
field can exist provided that £172 = —£271- One finds 
that a solution is given by 
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(3) 



When dealing with an interface separating a dielectric 
from a non lossy metal, e ^ 1 ^ 2 (2 yields a unique solu- 
tion to the problem. When accounting for losses in the 

material, , ei f 2 r is a complex number so that the dis- 
ci + £2 

persion relation cannot be solved using real K and real 
lu. It is necessary to consider a complex frequency or a 
complex wavector to find roots of the equation. Let us 
first choose 10 real. We denote K 8p the complex root of 
the equation ([3]). FigQJa) shows the dispersion curve 
obtained from the surface plasmon dispersion relation 
Eq. ([3]) when plotting to versus 9\eK sp . This curve ex- 
hibits a back-bending in the vicinity of the frequency of 
non-retarded surface plasmon u) p /y/2. The second possi- 
ble choice is to keep a real wavector K. We denote aj sp 
the complex root of the equation Eq. ©. Figfljb) shows 
the dispersion curve obtained when plotting d\tLu sp ver- 
sus K . It is seen in FigQJb) that this curve exhibits an 
asymptote for large wavevectors. 

Let us make two remarks regarding the dispersion re- 
lation. We first note that Eq.® is also a solution of 
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£i72 = e 27i which defines a zero of the reflection factor, 
i.e. the Brewster angle. It can be checked that the upper 
branch in FigQJb) is not a surface wave but the locus 
of the Brewster angle in the (u), K) plane. Finally, we 
note that a surface plasmon is a collective oscillation of 
charge density. When the frequency to is smaller than v e , 
the collective electron oscillation is ovcrdamped. In this 
frequency regime, Eq.Q describes a surface wave that 
has no longer the character of a surface plasmon. Thus, 
Eq.Q describes the Brewster angle for lo > ui p , a surface 
plasmon for lo p > lo > v ei a surface wave for v e > lo. 



III. GENERAL FIELD REPRESENTATIONS 

The aim of this section is to derive a general form of 
the surface plasmon field. To this aim, we first use the 
simple interface Green's tensor that yields the general 
form of the field for any given source distribution. We 
then extract the surface plasmon contribution which is 
defined to be the pole contribution to the Green's ten- 
sor. We will show that this procedure leads in a natural 
way to different representations that make use of either 
a complex wavevector or a complex frequency. We em- 
phasize that both representations will describe the same 
electromagnetic surface plasmon field E sp (a;, y, z, t). 

Let us suppose that an arbitrary source is located 
nearby the dielectric-metal interface. The electric field 
generated by the source j(r, t) is given by the relation 



E(r,i) = -fj, / dt' / dV G(r,r',t-t') 



^j(r',i') 



dt' 



(4) 



where /xo is vacuum permeability. A Fourier representa- 
tion can be written in the form 



G(r,r',t-t') = 



d 2 K f dio 



4tt 2 / 2tt 



• *g (K, z, z', lo) 

i[K(r-r')-w(t-i')] 



(5) 



Here, the integration variables lo and K Xl K y are real. 
The explicit form of the Green's tensor g(K.,z,z',Lo) in 
the presence of the interface is given in Appendix A. It 
is seen that the Fourier transform of the Green's ten- 
sor has poles given by the denominator of the Frcsncl 
factors for p-polarized field. For a dielectric/metal inter- 
face, they correspond to the surface plasmon as discussed 
previously. Thus, the Green's tensor can be split into 
two terms : the pole contribution that yields the sur- 
face plasmon and the remaining contribution that yields 
a regularized Green's tensor. 



G G re g -\- G S p, 



(6) 



where the pole contribution to the Green's tensor G sp can 

be explicitly derived using the residue theorem. G reg is 
the contribution of the regularized Green's dyadic. It can 
be shown that the Green's tensor can be evaluated us- 
ing a contour deformation in the complex plane and that 



the regularized term is essentially due to the contribution 
along the branch cut. This contribution is often termed 
cylindrical wave or creeping wave. The relative impor- 
tance of these terms is well documented in classical texts 
for radiowave a 30 i 31 . The analysis of their respective con- 
tribution was of practical importance in the early days 
of telecommunications as radiowaves were guided by the 
earth. This issue has been discussed recently in the con- 
text of optics^ 3 .. In this paper, we shall not pursue this 
discussion and focus instead on the surface wave contri- 
bution defined as the pole contribution. 
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dt' 
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G sp (r, r',t-t')- 



dt' 



(7) 

When solving Eq.Q, we can consider that lo is real 
and find a complex K sp or we can impose a real value to 
K and find a complex root lo sp . Thus, when extracting 
the poles, it is a matter of choice to consider that they are 
poles in the complex frequency plane or in the complex 
wavevector plane. We find cither a couple of poles lo sp 
and — uj* p or a complex wavevector pole K 2 p hence two 
poles for the component of the wavevector along the x 
and — -RTx.sp for a given component along the y 
It follows that we can cast 



axis -Kx,sp 
axis K v as K 2 



K 2 - 

x,sp 



K 2 

9 



<-y "sp 

the pole contribution to the Green's tensor in the form : 

%. p (K,z,z') J-„.JK,z,z') 
g sp (K,z,z ,lo) = — - — 1 , , — , (8) 



u — u) s 



where / Wsp (K, z, z') and /-^(K, z, z') are the residues 
of *g at w sp and respectively, or in the form : 

*g sp (K,z,z',uj) = 

fi<^ Bp (K y ,z,z',uj) f_ Kxap (K y ,z,z',uj) 



Ky^ ^-X, Sp 



+ 



where f K x , sp (K y ,z,z',u) and f^ Kx sp (K y ,z,z',u) are 
the residues of *g at if Xj sp and —K Xj sp respectively. 

These two choices leads to two different forms of the 
surface plasmon field given by Eq.([7]). We now examine 
these forms in detail. 



A. Surface plasmon field representation with a real 
wavevector 

In this section we derive the analytical form of the 
surface plasmon field using real wavevectors. For this 
purpose we evaluate the pole contribution to the Green's 
tensor by integrating in the complex lo plane. The com- 
plex pole w sp then yields a contribution for t — t' > that 
varies as exp(— iu> ap {t — t )). After integration, we find : 



G sp = H(t - t') 2^e 



d 2 K 



(2n) 2 
i[K-(r-r')-o) ip (t-t')] 



(-i)f w (K,z,z') 



(9) 
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where f u (H., z, z') is the residue of *g at uj sp . It is given 
in the appendix [A] It follows from Eq.([7|) that the field 
can be cast in the form of a linear superposition of modes 
with real wavevector and complex frequency : 



where K + 



2<Ke 



3 i(K-r+7 m |z|-w sp t) 



(10) 

where the amplitude i?(K,i) is given in appendix [Al 
n m = -z if z < and i if z > 0, and K = K/K. The 
surface plasmon field takes a form that looks as a mode 
superposition except that the amplitude i?(K,i) depends 
on the time t. Indeed, when describing a stationary field 
using modes that have an exponential decay, the ampli- 
tude is necessarily time dependent. In order to obtain a 
superposition of modes with fixed amplitudes, it is nec- 
essary to assume that all sources are extinguished after 
time t = so that we observe the field after it has been 
excited. In that case, the decay of the mode is well de- 
scribed by the imaginary part of lu sp . Eq. (| 1 0[) is thus 
well suited for fields excited by pulses. Note that the po- 
larization of each mode is specified by the complex vector 
K — y^n m , whose component along the z axis depends 
on the medium from which the field is evaluated. 



B. Surface plasmon field representation with a real 
frequency 

Let us now turn to the alternative choice. We con- 
sider the complex poles sp and —K x sp . The Green 
function can be cast in the form : 



G, 



2^ 



dK 
~2^ 

e iK *. sp(x-x') e iK y (y-y') e -iuj(t-t') 



~ fK Xtep (K y ,Z,z',U>) 



(11) 



if x — x' > 0, and : 

dto f dK. 



G sp = -i I - 
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f-K x , sp (K y ,z,z',uj) 



e -iK X: sp (x-x') e iK v (y-y') e -ioj(t-t') 



if 



< 



0. 



(12) 
and 



f if_ x , „ v (K y , z, z' ,ui) are the residues of 'g at A" Xi sp 
and — K Kt sp . They are given in appendix [X] When 
inserting this form in Eq.([7]), we again obtain a form 
for the field that is a superposition of modes whose 
amplitude depends on x : 



r duj 


fdK y 


1 2^ j 









E > (K y ,u;,x)(K+ -^n 



-n„ 



o«(--Kx, sp x+K y y+~/ m \z\-ujt) 



(13) 



(if x , sp x + K y y)/K sp and K 



(— .£sf XjSp x-|- K y y) /K sp . The amplitudes E > (K y , w, x) and 
E < (K y ,u>,x) are given in appendix lAl Again, it seems 
natural to have amplitudes of the modes that depend on 
x if one describes a homogeneous field using modes with a 
decay along x. A proper mode representation should use 
only fixed amplitudes. This is possible if all the sources 
lie in the x < region and the region of interest is the 
x > region. We then obtain a surface plasmon field 
that can be cast in the form : 



j'.(K-r+7 m |z|— u>t) 



(14) 



where K = JT X) sp x + K y y is complex and K = K/if sp . 
We conclude that stationary monochromatic fields with a 
finite size are well described by a representation that uses 
complex wavevectors and real frequencies. This equation 
is one of the main result of this paper. Indeed, it pro- 
vides a framework to develop surface plasmon Fourier 
optics. Similar representations have been postulated as 
ansatz to surface plasmons interferences^, propagation 
along a stripe^ or focussingi^. The framework introduced 
above provides a rigorous derivation of the form of the 
surface plasmon field valid in a region with no sources. 
Let us emphasize that this representation is well suited 
to discuss propagation for x > xo of a surface plasmon 
field known along a line x — xq. It is seen on Eq. (Tl4]) 
that propagation over a distance d amounts to multiply 
each mode by a factor exp(iK x d). In general, this in- 
volves modifying both the phase and the amplitude of 
the mode. Thus, it allows to discuss any surface wave 
diffraction problem. Finally, let us stress that this rep- 
resentation is valid for a complex wavevector K and a 
real frequency oj so that this representation is necessarily 
associated with a dispersion relation with backbending. 

To summarize, we have shown that the surface plasmon 
field can be represented using modes that have either a 
complex frequency or a complex wavevector. However, 
the amplitudes may still depend on either time or space. 
In the case of a field excited by a pulse, the represen- 
tation that uses a complex frequency is well suited. It 
is associated with the dispersion relation without back- 
bending. In the case of a stationary monochromatic ex- 
citation localized in space, a representation using modes 
with complex wavevectors is well suited. It corresponds 
to a dispersion relation with backbending. This simple 
analysis yields a simple prescription to choose the proper 
dispersion relation. Note that in the case of pulses lim- 
ited in space, both representations can be used. 



C. Surface plasmon field generated by a dipole 

For many applications, it is useful to know the field 
generated by a dipolc. For instance, when considering 
the field scattered by a subwavclcngth particle, the source 
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FIG. 2: A point-like dipole located nearby the dielectric- 
metal interface 



can be represented by an electric dipole. In addition, 
any source can be decomposed as linear superposition of 
dipolar sources. Here, we derive the surface plasmon field 
generated by a monochromatic point-like dipole charac- 
terized by its dipole moment po (see Fig. [2]). Note that 
there are other contributions to the field generated by the 
dipole at a distance typically smaller than a wavelength. 
The surface plasmon contribution is typically dominant 
for larger distance o 32 i 33 . 



For a vertical dipole p e 



-iuJot 



Poe 



; -iwotg located at a 



distance d below the interface, at x — y = 0, we obtain, 
using the cylindrical basis p, 6, z) : 



E, 



29k 



Hi 1} (K spP )p- 



; K sp rr(l) 



H^>(K sp p)n. 



M(K sp ,ujo) 



1 

eo 



.Ks 
Ti 



} (15) 



where and are Hankel functions of the first kind 
of zero-th and first order respectively, K ap is complex 
and verifies Eq. ^ with uj = u , M(K sp ,oJo) is given in 
appendix [B] and the other notations are defined above. 
The details of the calculation are given in appendix [B] 
Using the asymptotic forms of the Hankel functions, we 
obtain for the field of a vertical dipole, for p greater than 
a few 1/\K S J : 



E m = 29te 



iK sp p 



_^/K sp p 



K 

p- 7^f n m ) M' v {K ap ,w Q )pQ 



(16) 



where M' v {K sp , uj ) is given in appendix IBl Thus the sur- 
face plasmon field is analogous to a damped cylindrical 

iK sp p J£ 

wave ; __ with a polarization vector f — ^^-n m . For 

fir 



\J K sp p 

a dipole oriented along the x axis poe -1 " * 
we obtain : 



-iuiot-O- 



E, 
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gl (1) (^ p p) , 

K sp p 
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eo 



M(K sp ,u ) -poe^e^^e-*"*'} (17) 



where 9 = (x, p) . Using the asymptotic forms of the 
Hankel functions, we obtain : 
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M' h (K sp ,uj )po e 



(18) 



where M' h (K sp ,ujo) is given in appendix IB"1 The surface 
plasmon field is analogous to a damped cylindrical wave 

J , with the same polarization vector p — ^^-n m: 

V K spP lm 

times a factor cos 9, making the field vanishing in the 
direction perpendicular the dipole (here the y direction) 
and more intense in the dipole's direction (here the x 
direction). 



IV. DISCUSSION 

In this section, we revisit two fundamental issues in 
the field of surface plasmons : field confinment and large 
density of states. Confinment and supcrresolution are 
related to the existence of wavevectors with a modulus 
much larger than ui/c. In this respect the choice of the 
proper dispersion relation plays a key role as one has a 
cut-off wavevector whereas the other predicts no limit for 
the dispersion relation. Is there a limit to the resolution ? 
Is there a limit to the local density of states ? 



A. Super resolution 

Let us first discuss the issue of resolution when imag- 
ing with a surface plasmon driven at frequency u> by an 
external source. Recent experiments on far-field opti- 
cal microscopy^ launched a debat o 34 ' 35 about the role 
of surface plasmons in super-resolution imaging effects. 
In Ref. [2l| the dispersion curve with the asymptotic be- 
haviour has been invoked to stress the role of surface 
plasmons in the image formation with nano-resolution. 
The resolution was estimated to be X sp /2. Therefore, 
if the dispersion curve with the asymptotic behaviour is 
chosen, there seems to be no diffraction limit and only 
the amplitude decay of surface plasmon due to Ohmic 
losses in the metal limits the resolution. The effect 
of the back-bending of surface plasmon dispersion dis- 
cussed in Ref. [35[ limits the surface plasmon wavelength 
2n/9leK sp and therefore, the resolution. Clearly, both 
dispersion relations do not lead to the same conclusion 
and a prescription to choose one or the other is needed. 
Let us consider a situation where a surface plasmon is ex- 
cited locally by a stationary monochromatic field. From 
section 3, we know that it is valid to use a represen- 
tation with fixed amplitudes using modes with complex 
wavevectors and real frequencies. This implies that the 
dispersion relation with real frequency (with backbend- 
ing) is relevant. It follows that there is a cut-off spatial 
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frequency. Indeed, as K x may be complex, the prop- 
agation term exp(iK x x) introduces damping. In the 
case of a lossy medium, damping may be due to losses. 
However, even for a non-lossy medium (K sp is real), 
K x = (K 2 p — Ky) 1 / 2 can be imaginary. This occurs when 
K y exceeds the value K sp . This situation is the 2D analog 
of the evanescent waves with wavevector K larger than 
ui/c that cannot propagate in a vacuum. Clearly, K sp is 
a cut-off frequency and the propagation term exjp(iK x x) 
works as a low-pass filter that prevents propagation of 
fields associated with spatial frequencies larger than K sp . 
When dealing with lossy media, it is the real part of K sp 
that specifies the cut-off spatial frequency. It is seen in 
FigQ] that it is limited by the backbending of the disper- 
sion relation. 

In summary, when discussing imaging using stationary 
monochromatic surface plasmons, the relevant represen- 
tation is based on modes with a complex wavector and 
a real frequency given by Eq. (|14[) . This corresponds to 
a dispersion relation that has a backbending. It follows 
that the resolution is limited by the cut-off spatial fre- 
quency given by the maximum value of %KtK sp . 



B. Local density of states 

Let us now discuss the Local Density of States (LDOS). 
The density of states (DOS) is a quantity that plays a 
fundamental role in many domains. In particular, it al- 
lows to derive all thermodynamic properties of a system. 
In the case of an interface, the surface modes are con- 
fined close to the interface so that it is useful to introduce 
the Local Density of States (LDOS) that depends on the 
distance to the interfac o 36 i 37 . It allows to account for 
the huge increase of energy density close to an interface 
when surface waves are excite d 37 i 38 . It also plays a key 
role in defining the lifetime of a single emitter close to 
an interfac o 39 i 40 ' 41 i 42 ' 43 . In this context, the increase of 
the projected LDOS is usually normalized by the LDOS 
in a homogeneous medium (e.g. a vacuum) yielding the 
so-called Purcell factor. It is well-known in solid state 
physics that the density of states can be derived from the 
dispersion relation. More specifically, the DOS increases 
at a frequency uj when the dispersion relation is flat at 
that particular frequency. A quick look at FigfT] shows 
that different dispersion relations seem to predict differ- 
ent LDOS. While FigfjTb) predicts a very large peak at 
uj sp /\/2 due to the asymptote and no states above this 
frequency, Figflja) predicts a smaller peak and a non 
zero LDOS between u> sp /^/2 and u sp . Again, we see that 
a prescription is needed to choose the right dispersion 
relation. 

A standard procedure to derive the DOS in the recip- 
rocal space is based on the periodic boundary conditions. 
Assuming a surface of side L, the wavevector takes the 
form K = n x ^--k + riy2jjy. In the plane K x , K y , a mode 
has an area 4ir 2 /L 2 . It follows that the number of modes 



per unit area in cPK is given by cPK/47r 2 . When per- 
forming this analysis, both K x and K y are real. Thus the 
relevant representation uses real wavevectors and com- 
plex frequencies. The corresponding dispersion relation 
has no backbending and therefore presents a singularity. 
This is in agreement with another approach of the LDOS 
based on the use of the Green's tensor that predicts an 
asymptotic behaviour proportional to l/(z 3 |e + l| 2 ) 36 i 37 . 
Of course, this divergence is non physical. It is related 
to the modelling of the medium using a continuous de- 
scription of the metal. This model cannot be valid on an 
atomic scale. Before reaching the atomic scale, non-local 
effects must be taken into account. 



V. CONCLUSION 

The purpose of this work is to clarify several issues re- 
garding surface plasmons on flat surfaces. The first issue 
deals with the mode representation of the surface plas- 
mon field. We have shown that a surface plasmon field 
can be represented as a sum of modes with either a com- 
plex wavevector or a complex frequency. We have shown 
that a representation using complex frequencies is well 
adapted to fields excited by pulses and that a represen- 
tation using complex wavevectors is well adapted to sta- 
tionary monochromatic fields excited in a finite area. The 
latter representation provides a rigorous formula that can 
be used to analyse the diffraction of a stationary surface 
plasmon field. This should be very useful in order to de- 
velop a surface plasmon Fourier optics framework. This 
formula clearly shows that the maximum value of d\tK sp 
is a cut-off spatial frequency that gives an upper limit 
to the resolution or confinment that can be obtained us- 
ing surface plasmons. As a by product, we have derived 
the form of the surface plasmon excited by a dipole lo- 
cated below the interface. Finally, we have discussed 
how to choose the dispersion relation (with or without 
backbending) depending on the issue. To illustrate this 
procedure, we have shown that there is a resolution limit 
given by the maximum value of the wavevector at the 
backbending point. We have also shown that the local 
density of states should be analysed using the dispersion 
relation with a real wavevector. This yields a LDOS that 
diverges close to the interface in agreement with the re- 
sult obtained from the Green's tensor approach. 



APPENDIX A: CALCULATIONS OF THE 
SURFACE PLASMON FIELD 

For the plane interface system, it is convenient to use 
the representation due to Sipe 4 ^ that consists of a de- 
composition over elementary plane waves. We consider 
the interface such as the lower medium z < is denoted 
as medium 1 and the upper medium z > is medium 
2. We use the dyadic notation for the tensor. For in- 
stance, the s-componcnt of the electric field is given by 
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ssE = s(s • E). Using the notations of Eq. ([5]), we have 
for the field in the lower half-space z < and currents 
below z (z! < z) : 



*g(K, z, z', uj) = 5o(K, z, z', uj) 

+ o — [ §r ^ § + pr^pPi" 



-t-/lZ -l~f!Z 



(Al) 



where *go(K, z, z' , uj) denotes the Fourier transform of the 
Green's tensor of an infinite space filled by medium 1. 

For the field in the upper half-space (z > 0) and cur- 
rents still in the lower half-space (z! < 0), one has : 



*g(K, z, z', uj) = — [st s s + p 2 t p pt] e ~ 



271 Z 1^2 Z 



(A2) 



The Fresnel reflection and transmission factors are 
given by : 



7i ~ 72 
7i + 72 

27i 
7i + 72 



7162 - 72£i 
7i £ 2 + 72£i 

271^^62 

7i£2 +72£i 



(A3) 
(A4) 



7m = Y e m (^) — K 2 for m = 1,2 is chosen so that 
3m7, n > 0. This way, using e m (—u>*) = e m (w), one 
has y m (K*,—uj*) = — 7 m (if,o;). The square roots of 
dielectric constants y / e^T are chosen so that ^Kt^/e^l 0. 
s = K x z/K, pf = (Kz T j 1 K/K)/(k 0y /EI) and p 2 = 

(Kz - l2 K/K)/(k ^), fc = uj/c and K = Jk* + K$ 

is chosen so that 3mA > or 3mK = and EHeA' > 0. 
This way, one has p(— K, — oj*) = — p(K, u;)* when K is 
real and uj complex and p(— K*,— uj) = p(K, uj)* when 
K is complex and uj real. 

We extend the definition of g(K, z, z 1 , uj) to complex 
values of uj and assume that the denominator of the Fres- 
nel coefficients r p and t p (whose nullity is equivalent to 
Eq. ([3])) has two roots uj sp and -co* : 



C(K,uj) 



7i£2 + 72£i (uJ - uJ sp )(uJ + UJ* ) 



(A5) 



with 3muj sp < 0. *g then features two poles at uj sp and 
— uj* sp . The residues of *g at these poles can be calculated 

with z, z') = lim^^s [(uj — oj) g (K, z, z', w)] where 

cD denotes uj sp or -w* . It comes : 



f„jK,z,z') 



7i e 2 C(K,u! sp ) 
k 0y />T[ 2$ieuj sp 



A' 

(K n OT )p 1 " 

7m 



-1712 I7m|2| 



(A6) 



where n m denotes — z for z < and z for z > 0, 7 m 
denotes 7! for z < and 72 for z > 0. 71 and 72 depends 
on A and uj sp and pf on K and uj sp . 

Using e TO (-w*) = e* (w), 7 m (A",-w*) = -7 r * n (A,w), 
and pf (K, — cj*) = — pf (— K, w)*, we have : 



f- u: (K,z,z , ) = -f*(-K,z,z') 



(A7) 



Using Eq. J7J), ([9]) and (|A6[) . one can find the ampli- 
tudes in Eq. (TQ]) : 



£7(K,t) 



-/ip Ti^gl^v) /" d * T >e- iK *' 



dz'e- i7lZ 



diV^'pf. (A8) 



We now extend the definition of *g(K., z, z', uj) to com- 
plex values of K x and assume^ that the denomina- 
tor of the Fresnel coefficients r p and t p has two roots 

A x , sp = J Kip - K 2 and -A x , sp : 



7i £2 - 72£i 



7i£2+72ei e i - e l (Kk - A x , sp )(AT x + A x , sp ) 

(A9) 

with 3mA x , sp > 0, and where A sp depends on uj and 
is given by the dispersion relation. *g then features two 
poles at A Xi sp and — A Xj sp . 

The residues of *g at these poles can be calculated with 



fjr(K y ,z,z',u) = lim x ^ 



where A x denotes A x sp or — A X; sp 



(A x -A x )g(K,z,z', W ) 
It comes : 



« 7ie 2 7ie 2 - 7 2 ei 



fK x , p (K y ,z,z',uj) =- t 

ZK^ sp fcov £ i e i _ e 2 

(K+ - ^n m )p+(A x , sp x + K y y 7 uj)e~^ z 'e^ (A10) 



f-K X:Sp {K y ,z,z\uj) 



1 7ie 2 7i e 2-72£i 
2A XiS pA;oVeT e? - 4 



(K - ^n m )p+(-A XiSp x + A y y, 



ij!z' i-y m \z\ 



uj e 11 e 



(All) 



where K+ = *x,B P x + ir y y ^ = -A x , s x + A y y ; 

and the other notations are defined above. 

Using Eq. ([7|1. ([IT ]) . ([T2" ]) . (fATO)) and (TATT|) . one can 
find the amplitudes in Eq. (jT3|) : 



E > (K y ,uJ,x) = Mo 



1 71^2 7162 - 72£i 
2K X . sp fcoV^i e i ~ e i 



J dy'e-^y' J Az'e-^ 1Z ' J dt'e iut ' 



p+(A x , sp x + A y y, ^)-^(r', f) (A12) 



1 7i £ 2 7i £ 2 - 72£i 



E< (Ky,Uj, X )=„ - 2Kx spkoV _ 4 4 

/ dx'e- lK ^ x ' / dy'e- iK y y ' / dz'e" l7lZ ' / dt'e lwt ' 
p+ ( - A x ,s P x + K y y, u,). $L (r', t') (A13) 
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APPENDIX B: SURFACE PLASMON FIELD OF 
A DIPOLE 

The currents associated with the dipole poe~ luJat at 
a distance d above the interface are given by j(r, t) = 
2fHe [e- iuJot (-i)ujopo] S(r-(-d)z). Using the form of the 
surface plasmon field given by Eq. (|14|) . one can compute 
the amplitude with this expression and Eq. (|A12[) : 



E>{K v ,u) 



7i £ 2 - 72£i 7i £ 2 



p+. [2ttS(u> - w )p + 2ttS{uj + w )p*] (Bl) 



Hence, using Eq. JT4J) and the properties 
K(-K*) = -K*(K), K x , sp (K y ,-uj) = —K* sp (K y , ui), 
j m (K*,-u>) = -~(* n (K,u)), e m (-u>) = e* m (uj) and the 
definition of p+, it comes : 



E, 



-me 



i Jm \z\ -^pt Tifg -72£i 7i 72 

e e .2 Jl 



eo 



X : S P 



(B2) 



By derivating f^rj? = \H^\K sp p) with 

dK P ?K r - 

respect to x or y one finds f t \ y 4? K = 

Z7T i\ x ,sp 



iH[ l \K spP )p 

h£\k spP )]'pp 



A K iK r „ 

and KK 



r(i) 



J<spP 



2ir . 


^-x,sp 


p) ee 





Using these 



relations and the value of po, with M(K,co) = 
„ 7172 7i £2 - 72Ci ^ we get Eq iJ^ and (Tj^ According 

e l — e 2 

to Eq.(Hl]), the field found in Eq. (TS]) and CE7]) only 
apply in the x > half-space. Using symmetries 
arguments, one find easily that they also apply in the 
x < half-space. 

Eq. (fT5)) and (jTTJ) can then be simplified us- 
ing the asymptotical form of the Hankel functions 

H n (z) yjfe iZ -2™("+3). Denoting M' v (K sp , w ) = 

--y|e- l fM(X sp ,^ )^, we obtain Eq.©. Us- 
ing also the property of the Hankel functions 
H[(z) = H (z) - |i?i(z) and denoting M' h (K sp ,uj ) = 

-M(K sp ,ujo)-^, we obtain Eq. lfTB)) . 



2 -i* 
— p ^ 



ACKNOWLEDGMENTS 

This work was supported by the French National Re- 
search Agency (ANR) through Carnot Leti funding and 
by the French Ministry of Defense through a grant from 
the Direction generate de I'armement (DGA). 



* Electronic address: jean-jacques.grefTet@institutoptique.fr 

1 R.H. Ritchie, Phys. Rev. 106, 874 (1957). 

2 W. L. Barnes, A. Dereux, T.W. Ebbesen, Nature 424, 824 
(2003). 

3 T.W. Ebbesen, H.J. Lezec, H.F. Ghaemi, T. Thio, 
P.A. Wolff, Nature 391, 667 (1998). 

4 H. Liu, P. Lalanne, Nature 452, 06762 (2008). 

5 J.C. Weeber et al., Phys. Rev. B 64, 045411 (2001). 

6 S.I. Bozhevolnyi, V.S. Volkov, E. Devaux, J.Y. Laluet, 
T.W. Ebbesen, Nature 440, 04594 (2006). 

7 M. Quinten, A. Leitner, JR. Krenn, FR. Ausseneg, Opt. 
Lett. 23, 1331 (1998). 

8 Maier et al., Nature Mat. 2, 229 (2003). 

9 P. Andrew, W.L.Barnes, Phys. Rev. B 64, 125405 (2001). 

10 P. Anger, P. Bharadwaj and L. Novotny, Phys. Rev. Lett. 
96, 113002 (2006). 

11 S. Kuhn, U. Hakanson, L. Rogobete, V. Sandoghdar, Phys. 
Rev. Lett. 97, 017402 (2006). 

12 H. Rigneault, J. Capoulade, J. Dintinger, J. Wenger, 
N. Bonod, E. Popov, T.W. Ebbesen, P.F. Lenne, Phys. 
Rev. Lett. 95, 117401 (2005). 

13 J.B. Khurgin, G. Sun, R.A.Soref, Appl. Phys. Lett. 93, 
021120 (2008). 

14 C.E. Talley,J.B. Jackson, C. Oubre, NK Grady, CW Hol- 
lars, SM Lane, TR Huser, P. Nordlander, N.J. Halas, Nano 
Letters 5 1569 (2005) 

15 E. Prodan, C. Radloff, N.J. Halas, P. Nordlander, Science 
302, 419 (2003) 



K.L. Kelly, E. Coronado, L.L. Zhao and GC Schatz J. 
Physical Chemistry B 107 668 (2003) 
Z. Liu, J.M. Steele, W. Srituravanich, Y. Pikus, C. Sun, 
X. Zhang Nanoletters 5 1726 (2005) 

L. Feng, K.A. Tetz, B. Slutsky, V. Lomakin, Y. Fainman 
Appl.Phys.Lett. 91 081101 (2007) 

A.G Curto, F.J. Garcia de Abajo, Nano Lett. 8, 2479 
(2008) 

N. Fang, H. Lee, C. Sun, X. Zhang, Science, 308 534 (2005) 
I. Smolyaninov, J. Elliot, A.V. Zayats, C.C. Davis, Phys. 
Rev. Lett. 94, 057401 (2005). 

T. Taubner, D. Korobkin, Y. Urzhunov, G. Shvets, R. Hil- 
lenbrand, Science 313, 1595 (2006). 

T.V. Teperik, F. Garcia de Abajo, A.G. Borisov, M. Abdel- 
salam, P.N. Bartlett, Y. Sugawara, J.J. Baumberg, Nature 
Photonics 2, 299 (2008). 

J.J. Greffet, R.Carminati, K. Joulain, J.P.Mulet, S. Main- 
guy, Y. Chen, Nature 416, 61 (2002). 

F. Marquier, K. Joulain, J.P.Mulet, R. Carminati, J.J. Gr- 
effet, Y. Chen, Phys. Rev. B 69, 155412 (2004). 
M. Laroche, C. Arnold, F. Marquier, R. Carminati, J.J. 
Greffet, S. Collin, N. Bardou, J.L. Pelouard, Opt. Lett. 
30, 2623 (2005). 

R. Zia, M. L Brongersma, Nature Nanotechnology 2 426 
(2007) 

R.W. Alexander, G.S. Kovener, R.J. Bell, Phys. Rev. Lett. 
32, 154 (1974). 

ET. Arakawa, M.W. Williams, R.N.Hamm, R.H.Ritchie, 



9 



Phys. Rev. Lett. 31, 1127 (1973). 

30 A. Banos, Dipole Radiation in the presence of a conducting 
Half-space (Pergamon, NY, 1966) 

31 L. Felsen, N. Marcuwitz, Radiation and scattering of 
waves, (IEEE Press, NY, 1994) 

32 F. Pincemin, A. Sentenac, J.J. Greffet, J. Opt. Soc. Am. 
A 11, 1117 (1994). 

33 P. Lalanne, J.P. Hugonin, Nature Physics 2, 551 (2006). 

34 I.I. Smolyaninov, C.C. Davis, J. Elliot, A.V.Zayats, Phys. 
Rev. Lett. 98, 209704 (2007). 

35 A.Drezet, A.Hohenau, J.R. Krenn, Phys. Rev. Lett. 98, 
209703 (2007). 

36 K. Joulain, R. Carminati, J.P.Mulet, J.J. Greffet, Phys. 
Rev. B 68, 245405 (2003). 

37 K. Joulain, J.P. Mulet, F. Marquier, R. Carminati, J.J. Gr- 



effet, Surf. Sci. Rep. 57, 59 (2005). 

A.V.Shchegrov, K. Joulain, R. Carminati, J.J. Greffet, 
Phys. Rev. Lett. 85, 1548 (2000). 

R.R.Chance, A.Prock, R.Silbey, Adv. Chem. Phys. 37, 1 
(1978). 

K. Joulain, C. Henkel, Appl.Phys.B84 61 (2006) 
J. Nelayah, M. Kociak, O. Stephan, F. Garcia de Abajo, 
M. Tence, L. Henrard, D. Taverna, I. Pastoriza-Santos, L. 
M. Liz-Marzan, C. Colliex, Nature Physics 3, 348 (2007). 
F.J. Garcia de Abajo, M. Kociak, Phys. Rev. Lett. 100, 
106804 (2008). 

C. Chicanne et al, Phys. Rev. Lett. 88, 097402 (2002). 
J.E. Sipe, J. Opt. Soc. Am. B 4, 481 (1987). 
This is the case for a vacuum/metal (Drude model) inter- 
face. 



